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We present the study of weakly bound, neutron-rich nuclei using the nuclear shell model employing 
the complex Berggren ensemble representing the bound single-particle states, unbound Gamow 
states, and the non-resonant continuum. In the proposed Gamow Shell Model, the Hamiltonian 
consists of a one-body finite depth (Woods-Saxon) potential and a residual two-body interaction. We 
discuss the basic ingredients of the Gamow Shell Model. The formalism is illustrated by calculations 
involving several valence neutrons outside the double-magic core: ^-^°Ue and ^^-^^O. 

PACS numbers; 21.60.Cs, 21.10.-k, 24.10. Cn, 24.30. Gd 



I. INTRODUCTION 

The major theoretical challenge in the microscopic de- 
scription of weakly bound nuclei is the rigorous treatment 
of both the many-body correlations and the continuum 
of positive-energy states and decay channels. A fully 
symmetric description of the interplay between scattering 
states, resonances, and bound states in the many-body 
wave function requires a close interplay between methods 
of nuclear structure and nuclear reactions. This mutual 
cross-fertilization, which cannot be accomplished with- 
out overcoming a traditional separation between nuclear 
structure and nuclear reaction methods, is a splendid op- 
portunity for opening a new era in the nuclear theory of 
loosely bound systems. 

In many respects, weakly bound nuclei are much more 
difficult to treat theoretically than well-bound systems 
[1]. The major theoretical difficulty and challenge is the 
treatment of the particle continuum. For weakly bound 
nuclei (or for nuclear states above the particle threshold) , 
the continuum of positive-energy states and resulting de- 
cay channels must be taken into account explicitly. As 
a result, many cherished approaches of nuclear theory 
such as the conventional shell model (based on a single- 
particle basis of bound states) and the pairing theory 
must be modified. 

There are many factors which make the coupling to 
the particle continuum important. Firstly, even for a 
bound nucleus, there appears a virtual scattering into 
the phase space of unbound states. Although this pro- 
cess involves intermediate scattering states, the corre- 
lated bound states must be particle stable, i.e., they 
must have zero width. Secondly, the properties of un- 
bound states, i.e., above the particle (or cluster) thresh- 
old, directly reflect the continuum structure. In addition, 
continuum coupling directly affects the effective nucleon- 



nucleon interaction. 

The impact of the particle continuum was discussed 
in the early days of the multiconfigurational shell model 
(SM) in the middle of the last century. However, thanks 
to the tremendous success of the large-scale SM in terms 
of interacting nucleons assumed to be perfectly isolated 
from an external environment of scattering states [2~6] , 
the continuum-related matters had been swept under the 
rug. An example of an impact of the continuum that 
goes beyond the standard SM physics is the so-called 
Thomas-Ehrman shift [7,8] appearing in, e.g., the mirror 
nuclei ^^C, ^'^N, which is a salient effect of a coupling 
to the continuum depending on the position of the re- 
spective particle emission thresholds. The mathematical 
formulation of the problem of nuclear states embedded 
in the continuum of decay channels goes back to Fesh- 
bach [9], who introduced the two subspaces containing 
the discrete and scattering states. Feshbach succeeded in 
formulating a unified description of nuclear reactions for 
both direct processes in the short-time scale and com- 
pound nucleus processes in the long-time scale. As far 
as nuclear structure is concerned, the treatment of ex- 
cited states near or above the decay threshold has been a 
playground of the continuum shell model (CSM) [10-15]. 
Unfortunately, a unified description of nuclear structure 
and nuclear reaction aspects is much more complicated 
and became possible in realistic situations only at the 
end of the last century (see Ref. [16] for a recent review). 

In the CSM, including the recently developed Shell 
Model Embedded in the Continuum (SMEC) [17-19], the 
scattering states and bound states are treated on an equal 
footing. So far, most applications of the CSM, including 
SMEC, have been used to describe limiting situations 
in which there is coupling to one-nucleon decay chan- 
nels only. However, by allowing only one particle to be 
present in the continuum, it is impossible to apply the 
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CSM to 'Borromean systems' for which A- and (A-2)- 
nucleon systems are particle-stable but the intermediate 
(j4-l)-system is not. Various approaches, including the 
hyperspherical harmonic method or the coupled-channel 
approach, have been developed to study structure and re- 
ac;tion aspects of three-body weakly bound nuclei [20-22] . 
However, most of these models utilize the particle-core 
coupling which does not allow for the exact treatment of 
core excitations and the antisymmctrization between the 
core nucleons and the valence particles. 

The reason for limiting oneself to only one particle in 
the continuum in the CSM has been two-fold. First, the 
number of scattering states needed to properly describe 
the underlying dynamics can easily go beyond the limit of 
what present computers can handle. Second, treating the 
continuum-continuum coupling, which is always present 
when two or more particles are scattered to unbound lev- 
els, is difficult. There have been only a few attempts to 
treat the multi-particle case [23,24] and, unfortunately, 
the proposed numerical schemes, due to their complex- 
ity, have never been adopted in microscopic calculations 
involving multiconfiguration mixing. Consequently, an 
entirely different approach is called for. 

Recently, we formulated and tested the multiconfigu- 
rational shell model in the complete Bcrggren basis [25], 
the so-called Gamow Shell Model (GSM). (For applica- 
tion to two-particle resonant states, see also Ref. [26,27].) 
In this paper, GSM is applied to systems containing sev- 
eral valence neutrons. The single-particle (s.p.) basis of 
GSM is given by the Berggren ensemble, which contains 
Gamow states and the non-resonant continuum. The 
Gamow states [28] (sometimes called Siegert [29] or res- 
onant states) were introduced for the first time in 1928 
to study the a— resonances. Gamow defined complex- 
energy eigenstates E = Eq ~ ir/2 in order to describe 
the particle emission in the quasi-stationary formalism. 
Indeed, if one looks at the temporal part of such a state, 
which is e'-^"'/''e^'"*/'^^''\ one notices that the squared 
modulus of the wave function has the time-dependence 
oc e~^*, and one can identify fi/{r log 2) with the half-life 
of the system. 

Formally, the resonant states are generalized eigen- 
states of the time- independent Schrodinger equation with 
purely outgoing boundary conditions. They correspond 
to the poles of the ^-matrix in the complex energy plane 
lying on or below the positive real axis; they are regu- 
lar in origin and satisfy purely outgoing asymptotics. In 
the quasi-stationary approach with Gamow states, each 
observable O is complex. An interpretation of these com- 
plex values has been given by Berggren [30] : the real part 
of the matrix element gives the average value, while the 
imaginary part represents the uncertainty of the mean 
value. This is due to the finite lifetime of the Gamow 
state which implies that none of the measurements in 
this state can have a well-defined probability. 

In the previous pilot work [25] we showed first appli- 
cations of the GSM. In this work, we give the details of 
calculations and demonstrate first applications of GSM 



to particle distributions and transition matrix elements. 

The paper is organized as follows. Section II discusses 
how to calculate the matrix elements in the Berggren 
basis. The completeness relations valid for the single- 
particle resonant states are briefly reviewed in Sec. Ill, 
and some numerical examples involving the Berggren set 
of the Woods-Saxon potential are presented. Section IV 
describes the GSM Hamiltonian used in our work. The 
extension of the completeness relations to the many-body 
case is described in Sec. V. Sections VI and VII contain 
the GCM analysis of is-^sq and ^-'^°Re, respectively. 
Finally, Sec. VIII contains the main conclusions of the 
paper. 

II. MATRIX ELEMENTS IN THE BERGGREN 
BASIS 

Gamow functions are solutions of the Schrodinger 
equation which arc regular at the origin and have the 
outgoing wave asymptotics; i.e., the radial part behaves 
as e'*"' at large distances. In the case of the spherical 
one-body potential, the resonant wave function (f)nji car- 
rying the s.p. angular momentum £j can be written as 
a product of the usual angular part and the radial wave 
function Unje{r)/r. It is customary to introduce the no- 
tation 

Unjiir) = u„ji{r)* (1) 

4'njl = 4>njl>{u ^ u). (2) 

For bound states, one can always introduce a phase con- 
vention which makes the radial wave function real. That 
is, for bound states = (/)„j^. The following discus- 
sion concerns the specific properties of the radial Gamow 
wave functions Unje{r). Of course, one should always re- 
member that the angular part is always present, but its 
treatment is standard. Consequently, in the GSM, only 
radial matrix elements require special attention. 

A. Normalization of Gamow states and one-body 
matrix elements 

The norm Ni of a resonant state, 

Nf = / u^{r)dr, (3) 
Jo 

and the radial matrix elements calculated in the Berggren 
basis, 

r+oo 

Oif= Uf{r)0{r)ui{r)dr, (4) 
Jo 

are diverging, but this difficulty can be avoided by means 
of a regularization procedure [31-36]. Zel'dovich pro- 
posed to multiply the integrand of a radial matrix ele- 
ment by a Gaussian convergence factor [31]: 
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{uf\0\ui) 



lim 



+00 

e" 



Uf{r) 0{r) Ui{r) dr. (5) 



In this expression, \uf) and stand for single-particle 
(s.p.) states, and 0(r) is a radial part of a one-body 
operator. Using the Zcl'dovich regularization method, 
Berggren has shown [34] that the Gamow states, together 
with scattering states, form a complete basis. In partic- 
ular, using definition (5) with 0(r)=l, one can demon- 
strate that all Gamow states can be orthonormalized, 
and such orthonormalized functions can be used to cal- 
culate matrix elements. Unfortunately, the method of 
Zel'dovich, even though important on formal grounds, 
cannot be used in numerical applications due to the dif- 
ficulty in approaching the limit in (5) for diverging inte- 
grals. 

An equivalent and more practical procedure, justified 
by the apparatus of the analytical continuation, was pro- 
posed by Gyarmati and Vertse [35]. For that, let us define 
the following functional on (—00; Vum) [37]: 



FiVo) 



f 



(6) 



where Vo is the depth of the potential generating s.p. 
wave functions Uf and Ui, Vum is the depth of the po- 
tential for which one of these functions is bound and the 
other one is at zero energy, and 0{r) is some analytical 
operator. This functional is defined in such a way because 
the integral converges in the domain (—00; Vum). It rep- 
resents the radial matrix clement (w/|0|ui)/||u/||/||ui|| 
between two not necessarily normalized wave functions. 
Since Uf,Ui are bound, one can make them real, = Uf, 
u* = Ui. In Rcf. [37], the analytical continuation of F 
is made using the Padc approximants. In the present 
work we shall refer to the technique of the complex ro- 
tation [35] which allows calculation of F with Vo > Vum- 
To see that, let us call /(r) one of three integrands 
Uf{r)0{r)ui{r), u'j{r), orwf(r) and let us take Vo < Vum- 
Since / is analytical on C (see Fig. 1), then, following 
the Cauchy theorem, one has: 



fiz)dz+ / f{z)dz 



+ 



[ f{z) 



dz = 0. 



(7) 



Since / decreases exponentially for Re[z] > 0, the in- 
tegral j^^ f{z) dz ^ if Rf ^ +00. For the same rea- 
son, the integrals J^^ f{z) dz and J^^ f{z) dz converge if 
Rf +00. Consequently, for Rf — > -|-oo one obtains: 



+00 



f{r) dr 



r+00 
Jo 



X ■ e'^Je'" dx. 



(8) 



Hence, on the interval {—(x;Viim), one can define F 
by Eq. (6), with the norm (3) given by 



Im[rJ 




R+Rj 



Re[r] 



FIG. 1. The path in the complex coordinate space corre- 
sponding to the complex rotation by angle 9. R is the point 
from which the exterior complex rotation starts. R is large as 
compared to the nuclear radius; hence it is assumed that the 
nuclear potential is negligible for r>R. 




u1{r) dr 



+00 



uf{R + x-e'^)e'^dx (9) 



and with the matrix element (4) of the form: 



Uf{r) 0{r) Ui{r) dr 



+ 



/ [uf{R + x-e'^)0{R + x-e'^) 
Jo 

xui{R + x-e''^)e''' dx] . 



(10) 



li Uf and Ui arc bound- or decaying-state wave func- 
tions, one can write fc/ = |fc/|e~*"^ and ki = |fci|e^*"'. 
As Uf{z) ~ af{z)e^^''' and Ui ^ ai{z)e^''^^ when Re]^;] — > 
-hoc, with 0/ and Oj the algebraic increasing functions, 
the integrals defining F converge if one takes 



6 > aj + ai. 



(11) 



In addition, the expression for F is analytical because F 
is a function of converging integrals of analytical func- 
tions. Square roots in Eq. (9) cause no problems because 
Nf and have always a positive real part. Conse- 
quently, following the theorem of analytic continuation, 
Eq. (6) defines also F foY V > Vum . In this way, one may 
calculate the radial matrix elements of resonance states 
which a priori are not normalizable. 



B. Scattering states 

Scattering states represent the non-resonant contin- 
uum and explicitly enter the completeness relations dis- 
cussed in Sec III. Their asymptotic behavior at r — > -|-oo 
is: 

u{r) ~ C+H^+\kr) + C-H^,-\kr), (12) 
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where Hjf^^ denote Hankel (or Coulomb) functions. As 
usual, u(r) is normalized to the Dirac (5— distribution: 



/. + 00 

/ U{ 

Jo 



u{k, r) u{k', r) dr = d{k — k'), 



which gives 



c+c- = 



2tt' 



(13) 



(14) 



Knowing u{r), , H^^\ and their derivatives at point 
R, one may determine coefficients C+ and C~ up to a 
normalization factor by solving the set of linear equa- 
tions: 



u{R) = C+H^,+\kR) + C-H^,-\kR) 



(15) 



u'{R) = kC+ 



(+) 
e,ri 



dz 



+ kC- 



z=kR 



dH. 



(-) 



dz 



z=kR 



Finally, the scattering state is normalized to satisfy the 
condition (14). 



C. Matrix elements involving scattering states 

The calculation of matrix elements involving the scat- 
tering states is based on the complex rotation (9,10). 
However, the analytical continuation should be intro- 
duced differently than for resonant states. The one-body 
matrix element can be written as: 



F{kf) 



Ufir)V{r)u,{r) dr + AfA,F++{kf) 

AfB,F+-{kf)+BfAiF-+{kf) (16) 
BfB,F__ikf), 



where 



Uf = Afu'^ + BfUj ; 



• Ui = AiUf 



BiU~ , where hi in F is fixed and, 
in general, Ui can be either a bound, resonant or 
scattering state; 

r + oo 

• FsfsAkf)= / uyiR + x)0{R + x)Ui'{R + x)dx 
with s/, Si G (-I-, -). 

This separation is necessary because the presence of in- 
coming and outgoing waves in the same integral does not 
allow one to find a unique path in the complex plane 
along which the integrand decreases exponentially. Con- 
sequently, for each Fs^si one has to consider the domain 
of the complex plane where it converges, and then one 
performs an analytical continuation with the appropri- 
ate angle 9sfSi- 



Certain integrals cannot be regularized in the above 

sense. Those include F^ and F ^ with Ui = uj. For 

0(r) = 1, the integrand tends toward a constant value 
at -|-oo, independently of the value 9+-- This can be im- 
mediately seen for neutrons, because with z = R + x- e^^ 
and with \z\ +oo, the product u^{z) ■ u~{z) — > 
const X e^''^ x e~'^'^^ = const, and the corresponding inte- 
gral diverges. In this case, however, it is easy to see that 
the integral is in fact a 5— distribution, and it can be 
calculated by using a discrete representation of the Dirac 
5— function. 



5{k - ko) 



Ok,ko 



(17) 



with Afc being the discretization step in k. 



III. COMPLETENESS RELATION INVOLVING 
SINGLE-PARTICLE GAMOW STATES 

There exist several completeness relations involving 
resonant states. As shown by Lind [38], they all can 
be derived from Mittag-Lcffler theory. In the following, 
we briefly discuss the Berggren completeness relation [34] 
which is used in our paper. The following discussion will 
concern the s.p. radial wave functions corresponding to 
a given partial wave 

We begin from the completeness relation of Newton 
[39]: 



r+oo 

(wn I + / \uk) {uk I dk 
„ Jo 



(18) 



where |u„) are the normalized bound states and \uk) 
are the scattering states along the real energy axis nor- 
malized according to (13). In the basis (18), one can 
expand any bound state or scattering state with real 
energy. Unfortunately, in the presence of narrow reso- 
nances, the discretization of the real energy continuum 
becomes cumbersome. The complex-energy formalism of 
Gamow states offers a simple remedy to this difficulty. 

In order to derive the completeness relation with 
Gamow states, one has to deform the integration contour 
into the complex fc— plane, as shown in Fig. 2. Following 
the residuum theorem, one obtains: 

- j \uk){uk\dk -\- I \uk)f^\dk 
Jo Jl+ 

= 2in^Res{\uk){^k\)k=k^, (19) 

k„ 

where fc„ are the poles of \uk){uk\ lying between the real 
axis and the complex contour. 

In general, the scattering wave function Uk can be writ- 
ten as : 
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FIG. 2. Representation of the complex fc— plane, showing 
the positions of bound states, resonances, antiresonances, and 
the antibound states. L+ is the contour representing the 
non-resonant continuum. The Berggren completeness rela- 
tion involves the bound states, decaying states lying between 
L+ and real fc— axis, and the scattering states on L+. The 
contour L+ has to be chosen in such a way that all the poles 
in the discrete sum in Eq. (28) are contained in the domain 
between L+ and the real energy axis. 



In the above expression, stands for the Jost function 
[40]: 



. du du 
dr dr ' 



(21) 



-oo, and u 



where ^ H^^{kr) when r - 
Au^ + Bu^ . For bound and decaying states, J^{k)^Q 
{J~{k) = only for capturing states with the incoming 
wave asymptotics) . Consequently, the pole of Uk{r) cor- 
responds to zero of J'^{k). When fc ^ fc„, i.e., when U]^ 
approaches the resonant state, then: 



\ / ^1 (fcn) I + I / +1 

The derivative of the Jost function dX k — k„ is 



d£ 

dk 



ij {kn)Reg 



k=k„ 



(r) dr 



As this derivative is nonzero, we have as fc — > A:„: 

'dj+' 



J+{k)^{k-K,) 



dk 



k=k„ 



hence 



\uk){uk\ 



1 



27ri(fc — fc„) 
where the normalized resonant state is 



(22) 



(23) 



(24) 



(25) 



Un{r) = ul (r) Reg 



+ 00 



uj^ (r) dr 



Finally, the residuum at fc = kn is 

1 



Res{\uk){uk\)k=kr, 



2iTT 



Un){Ur. 



and the completeness relation follows immediately: 



^ \Un){Ur, 



\uk){uk\dk = 1. 



(26) 



(27) 



(28) 



In the above equation |u„) are the Gamow states (both 
bound states and the decaying resonant states between 
the real fc— axis and the complex contour). Relation (28) 
is the Berggren completeness relation which allows one to 
expand the states with complex fc inside the zone between 
real fc— axis and the complex contour. One may notice 
again that the resonances in Eq. (26) are normalized us- 
ing the squared wave function and not the modulus of 
the squared wave function. This is a consequence of the 
analytical continuation which is used to introduce the 
normalization of Gamow states. 

Figure 2 illustrates the ingredients entering Eq. (28). 
The resonant states, the poles of the S matrix, are repre- 
sented by the dots. They are divided into the bound, 
decaying, capturing, and antibound states (see, e.g., 
Refs. [34,41,38]). The relation (28) involves the bound 
and decaying states and the contour _L+ lying in the 
fourth quadrant of the complex-fc plane. 

In practical applications, one has to discretize the in- 
tegral in (28) [42,43]: 



\uk){uk\dk ~ ^ \ui){ui\, 



(29) 



where Ui{r) = -y/Afe.Ufc. (r) and Afc^ is the discretization 
step. It follows from the definition of Ui{r) that 

{ut\uj) = Si^j, (30) 

and the discretized Berggren relation (28) takes the form: 



\Un){Un \ + ^ \Ui){Ui 



(31) 



This relation is formally identical to the standard com- 
pleteness relation in a discrete basis and, in the same way, 
leads to the eigenvalue problem — E\'^). However, 

as the formalism of Gamow states is non-hermitian, the 
matrix H is complex symmetric. 

Up to this point, the choice of the contour in Eq. (2) 
has been completely arbitrary. In practice, however, one 
wants to minimize the number of discretization points 
Nd along L+. This can be achieved if the scattering 
functions on the contour (or, rather, their phase shifts) 
change smoothly from point to point. This condition 
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can be met if the contour does not lie in the vicinity 
of a pole, especially the narrow resonant state. If this 
condition is met, the states appearing in the basis (31) 
can be naturally divided into: 

• Bound states ~ lying on the imaginary fc-axis (or 

negative real energy axis), 

• Narrow decaying states - lying close to and below 
the real fc-axis (or below the positive real energy 
axis). Those states can be interpreted as physical 

resonances of the system, 

• Non-resonant continuum - represented by the scat- 
tering states along Physically, those are the 
building blocks of the non-resonant background. 

These definitions can be extended to many-body states 
in the complex energy (or momentum) plane. In the 
following, we shall clearly distinguish between resonant, 
resonance, and non-resonant states. 

The completeness relations derived above hold in every 
channel. Consequently, in practical calculations, 
one has to take different contours for different partial 
waves. As discussed below, the choice of the contour de- 
pends on the distribution of resonant states in the com- 
plex fc- plane. 

In a number of papers (see, e.g., Refs. [44-46]), reso- 
nant states were applied to problems involving continuum 
in the so-called pole expansion, neglecting the contour in- 
tegral in Eq. (28). The importance of the contour con- 
tribution was investigated in Refs. [47,48,25-27] where it 
was concluded that if one is aiming at a detailed descrip- 
tion, the non-resonant contribution must be accounted 
for. This point will be clearly seen in several examples 
discussed below. 



^3 = 2.0 — iO.O. The contour is discretizcd with a differ- 
ent number of points: n=60, 80, 100, 120, 160, and final 
results are obtained using the Richardson extrapolation 
method. In the examples considered in this section, we 
shall expand the 2^3/2 state, |wws)) either weakly bound 
or resonant, in the basis |u'wsB(fc)) generated by the WS 
potential of a different depth: 

|mWs) = y^CfcjMwSB(fci)) 

i 

+ / c{k)\uvf SB {k))dk, (33) 

cf. Eq. (28). In the above equation, the first term in 
the expansion represents contributions from the resonant 
states while the second term is the non-resonant contin- 
uum contribution. Since the basis is properly normalized, 
the expansion amplitudes meet the condition: 

[ c\k)dk = l. (34) 

In all cases considered, the OP3/2 and IP3/2 orbitals are 
well bound (by ~ 40 MeV and ~ 18 MeV, respectively) 
and do not play any role in the expansion studied. 

In the first example, we shall expand the 2p3/2 s.p. res- 
onance (0.25-10.20 MeV) of a WS potential of the depth 
Vo=58 MeV in the basis generated by the WS potential 
of the depth Vg^=59 MeV (here the 2^3/2 s.p. resonance 
has an energy of 0.16- iO. 09 MeV). The density of the ex- 
pansion amplitudes is shown in Fig. 3. One can see that 
the contribution from the non-resonant continuum is es- 
sential even though the 2p3/2 basis state is a resonance. 
In this example, as one might expect, the contribution 
from the resonance state in the basis is dominant. 



A. Completeness of the one-body Berggren basis: 
illustrative examples 

In this section, we shall discuss examples of the 

Berggren completeness relation in the one-body case. 
The s.p. basis is generated by the spherical Woods-Saxon 
(WS) potential: 



V{r) = -Vof{r)-V,o^l-s-'^^^''^ 



fir) 



1 -|- exp 



r-Ro 



r dr ' 
1 



(32a) 
(32b) 



In all examples of this section, the WS potential has the 
radius Ro = 5.3 fm, diffuseness d = 0.65 fm, and the spin- 
orbit strength Vso = 5.0 MeV. The depth of the central 
part is varied to simulate different situations. 

The complex contour corresponds to three straight 
segments in the complex fc— plane, joining the points: 
ko = 0.0 - iO.O, ki = 0.2 - iO.2, fca = 0.5 - iO.O and 




0.2 0.3 
Re(k) 

FIG. 3. Distribution of the squared amplitudes c^{k) of 
the s.p. state 2^3/2 of one WS potential (Vo = 59 MeV) in 
the s.p. basis generated by another WS potential (Vq =58 
MeV). The amplitudes of both real (solid line) and imagi- 
nary (dashed line) parts of the wave function are plotted as a 
function of Rc[fc]. The height of the arrow gives the squared 
amplitude of the 2p3/2 resonance contained in the basis. 



The second example shown in Fig. 4 deals with the 
case of a 2^3/2 state that is bound in both potentials. 
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Here Vo=64 MeV and Vf\ '=62 McV, and the 2^3/3 state 
lies at -0.91 MeV and -0.33 MeV, respectively. As in the 



■ Vp=64 MeV, 2p3/2 : -0.91 MeV ■ 

0.6 : Vg=62 MeV, 2p3,2 : -0.33 MeV : 




real part 
imaginary part 



-0.2 0.0 0.2 0.4 0.6 0.8 
Re(k) 

FIG. 4. Similar as in Fig. 3 but for the bound 2p3/2 s.p. 
state of tlie WS potential with V(i=64 MeV expended in the 
basis generated by another WS potential {Vq^^ = 62 MeV). 
The height of the arrow gives the squared amplitude of the 
bound 2p3/2 state at the value of -Im[fc] (the corresponding 
A;— value is purely imaginary). 

previous case, the contribution from the resonant (here: 
bound) state in the Berggren basis dominates and the 
contribution of the non-resonant continuum is small al- 
though non-negligible. In this figure, one can notice the 
small cusp at Re[fc]=0.2, even though the density c^{k) is 
an analytic function of k. This apparent paradox is due 
to the fact that c^(fc) is plotted as a function of Rc[fc] 
(Im[A;]=0). Moreover, the path in the complex plane is 
continuous but not derivable at k = 0.2 — iO. 2. These two 
aspects contribute to the appearance of the 'discontinu- 
ous feature', which of course has no physical meaning. 

An interesting situation is presented in Fig. 5. Here the 
unbound 2p3/2 state (Vo=59 McV) is expanded in a WS 

basis containing the bound 2p3/2 level {Vq —62 MeV). 
Consequently, the non-resonant continuum has to supply 

the imaginary part of the resonance's wave function. The 
last example (Fig. 6) corresponds to Vo=64 MeV and 

Vq —59 MeV. This is the most intriguing case since one 
expresses a bound (real) state in the basis which contains 
only complex wave functions (the contribution from well 
bound 0p3/2 and 1^3/2 s.p. states is negligible). In this 
case, the non-resonant continuum annihilates the imag- 
inary component of the 2p3/2 s.p. resonance contained 
in the basis. Indeed, in this example, the contribution 
from the contour is dominant. To see the convergence of 
the wave function obtained by the expansion method, in 
Fig. 7 we show the root mean square (RMS) deviation of 
the calculated 2^33/2 wave function from the exact result 
as a function of the number of discretization points along 
the contour. While the wave functions converge fairly 
quickly, the convergence of complex energies is slightly 
slower; hence, one has to employ the Richardson extrap- 
olation method to get an energy precision of the order of 
a keV. 



V„=59MeV, 2p3,2:0.16-i0.09McV : 
Vo=62 MeV, 2p3,2 : -0.33 MeV 
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FIG. 5. Similar as in Fig. 3 except for the 2^3/2 resonance 
of the WS potential with Vb=59 MeV expended in the basis 
generated by another WS potential (VJ, = 62 MeV). 



IV. GAMOW SHELL MODEL HAMILTONIAN 

The GSM Hamiltonian applied in this work consists 
of a one-body term and a zero-range two-body interac- 
tion. The spherical one-body potential was taken in a WS 
form (32a). In our study, resonant states are determined 
using the generalized shooting method for bound states 
which requires an exterior complex scaling. The numer- 
ical algorithm for finding Gamow states for any finite- 
depth potential has been tested on the example of the 
Poschl-Teller-Ginocchio (PTG) potential [49], for which 
the resonance energies and wave functions are known an- 
alytically. Energies of all PTG resonances with a width 
of up to 90 MeV are reproduced with a precision of at 
least 10"^ MeV. 

Contrary to the traditional shell model, the effective 
interaction of CSM cannot be represented as a single 
matrix calculated for all nuclei in a given region. The 
GSM Hamiltonian contains a real effective two-body 
force expressed in terms of space, spin, and isospin coor- 
dinates. The matrix elements involving continuum states 
are strongly system-dependent, and they have to be de- 
termined for each case separately. This creates an ad- 
ditional difficulty, but there is also a pay-off. Namely, 
the resulting two-body matrix elements fully take into 
account the spatial extension of s.p. wave functions. 

In this work, as a residual interaction we took the sur- 
face delta interaction [50] 

V{1, 2) = -VsmS {n - r*2) ^(ri - R) (35) 

with the the same value of R as in the WS potential. In 
our exploratory GSM calculations, we consider two cases: 
(i) the chain of oxygen isotopes with the inert ^^O core 
and active neutrons in the sd shell, and (ii) the helium 
chain with the inert *He core and active neutrons in the 
p shell. The parameters of the GSM Hamiltonian are 
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Vg=64 MeV, 2p3/2 : -0.91 MeV 
Vo=59 MeV, 2py^ : 0.16-i0.09 MeV 




real part 
■ imaginary part 



0.2 0.4 0.6 0.8 1.0 

Re(k) 

FIG. 6. Similar as in Fig. 3 except for the 2p3/2 bound 
state of the WS potential with V(3=64 MeV expended in the 
basis generated by another WS potential (Vj}^^=59 MeV) 
containing a resonance. 



Variant 


R (fm) 


d (fm) 


Vo (MeV) 


v;o (McV) 


VsDI 




(frn) 


(fm) 


(McV) 


(McV) 


(MeVfm^) 




3.05 


0.65 


55.8 


6.06 


700 


'Be 


2.00 


0.65 


47.0 


7.50 


1670 



TABLE I. Parameters of the GSM Hamiltonian used in the 
calculations for the oxygen isotopes ("^'"O" parameter set) 
and the helium isotopes ("^Hc" parameter set): WS radius 
R, WS diffuseness d, WS strength Vo, spin-orbit strength V^so, 
and strength of the residual SDI interaction Vsdi- 



summarized in Table I. The WS potentials have been 
adjusted to s.p. states in one-neutron nuclei ^''O and 
^He. For ^'^O, the resulting WS 0^5/2 and lsi/2 states 
are bound, with s.p. energies -4. 142 MeV and -3.272 MeV, 
respectively, and 0^3/2 is a resonance with the s.p. en- 
ergy 0.898-10.485 MeV. The agreement with experimen- 
tal data (e^^P+ =-4.143 MeV, e^^^+ = -3.273 MeV, and 

e^^^ = 0.942 MeV, = 96 keV) is excellent. The 

strength of the SDI has been adjusted for a given con- 
figuration space to the experimental two-neutron sepa- 
ration energy of "'^^O. Since only 0^3/2 is a s.p. reso- 
nance, we shall include only a ^3/2 non-resonant contin- 
uum. In fact, the completeness relation requires taking 
the non-resonant continua corresponding to all partial 
waves {£, j). However, if for a given partial wave no reso- 
nances arc included in the basis, the corresponding non- 
resonant continua can be chosen along the real momen- 
tum axis. Since, to the first order, the inclusion of these 
continua should only result in the rcnormalization of the 
effective interaction, they can be ignored in most cases, 
except for Sec. VB. 

The nucleus ^He, with one neutron in the p shell, is 
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Number of non-resonant continuum points 
FIG. 7. The root mean square (RMS) deviation of the 2p3/2 
wave function of the WS potential (Vo = 64 MeV) obtained by 
a diagonalization in a basis generated by another WS poten- 
tial {V(f=59 MeV) from the exact wave function (obtained by 
a direct integration of the Schrodinger equation). The RMS 
deviation is shown as a function of the number of discretiza- 
tion points along the 2p3/2 contour. 



unstable with respect to the neutron emission. Indeed, 
the J'' = 3/2^ ground state of ^He Ues 890 keV above 
the neutron emission threshold and its neutron width 
is large, r=600kcV. The first excited state, 1/2 is 
a very broad resonance (r=4MeV) that lies 4.89 MeV 
above the threshold. Our WS potential yields single- 
neutron resonances and pi/2 at £^=0.745-10.32 MeV 
and ii;=2. 130-12.936 MeV, respectively. In our model 
space we take resonances 0^3/2, Opi/2, and the two as- 
sociated complex continua p^/2 and pi/2- The strength 
of the SDI has been adjusted for a given configuration 
space to the experimental two-neutron separation energy 
of ^He. 

For the A'^-body problem, the Hamiltonian matrix con- 
tains one- and two-body matrix elements. The one-body 
part corresponds to s.p. energies of basis states and con- 
tributes only to diagonal matrix elements. In general, 
the calculation of two-body matrix elements is performed 
by splitting the radial integral into 16 terms correspond- 
ing to all different possible asymptotic conditions of s.p. 
wave functions. Then, each term is regularized separately 
by an appropriate choice of angle of the external complex 
scaling, cf. Sec. II. 



V. MANY-BODY COMPLETENESS RELATION 
WITH GAMOW STATES 

The discretizcd basis (31) can be a starting point for 
establishing the completeness relation in the many-body 
case, in a full analogy with the standard shell-model in 
a complete discrete basis, e.g., the harmonic oscillator 
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basis. In this case one has: 

5^i*„)($;i^i, 



(36) 



where the A^-body Slater determinants l^fn) have the 
form \(f)i ...0jv), where \(j)k) are resonance (bound and 
decaying) and scattering (contour) s.p. states. The ap- 
proximate equality in (36) is an obvious consequence of 
the continuum discretization, similarly as in (31). Like 
in the case of s.p. Gamow states, the normalization of 
the Gamow vectors in the configuration space is given by 
the squares of shell- model amplitudes: 



(37) 



and not to the squares of their absolute values. 

In the particular case of two-particle states, the com- 
pleteness relation reads: 



(38) 



This relation can be used to calculate the two-body ma- 
trix elements. 



A. Determination of many-body bound and 
resonance states 

Before discussing completeness relations in the many- 
body case, let us describe the method of selecting many- 
body resonances. In a standard shell model, one often 
uses the Lanczos method to find the low-energy eigen- 
statcs (bound states) in very large configuration spaces. 
This popular method is unfortunately useless for the de- 
termination of many-body resonances because of a huge 
number (continuum) of surrounding many-body scatter- 
ing states, many of them having lower energy than the 
resonances. A practical solution to this problem is the 
two-step procedure proposed in Ref. [25]: 

• In the first step, one performs the pole approxi- 
mation, i.e., the Hamiltonian is diagonalized in a 
smaller basis consisting of s.p. resonant states only. 
Here, some variant of the Lanczos method can be 
applied. The diagonalization yields the first-order 
approximation to many-body resonances 
where index i {i = 1, . . . ,N) enumerates all eigen- 
vectors in the restricted space. These eigenvectors 
serve as starting vectors (pivots) for the second step 
of the procedure. 

• In the second step, one includes couplings to non- 
resonant continuum states in the Lanczos subspace 
generated by {j e [1, . . . ,N]). 

• Finally, one searches among the M solutions \^j:k), 
{k = 1, . . . , M) for the eigenvector which has the 
largest overlap with 



This procedure is a variant of the Davidson method. Ob- 
viously, in the search for bound states, both Lanczos and 
Davidson methods can be used. In the following, we shall 
show that the procedure outlined above allows for an ef- 
ficient determination of physical states within the set of 
all eigenvectors of a given Lanczos subspace. 

As a representative example, let us consider the cases 
of and ^^O with the core of ^^O, i.e., two- and 
four-particle systems, respectively. Here wc employ the 
complex (i3/2 contour corresponding to two segments de- 
fined by the points: fci = -|- iO, = 0.2 — i0.05, and 
^3 = 0.4 — iO.O. The contour is discretized with 17 and 
9 points for and ^"O, respectively. From all the pos- 
sible many-body configurations, we only keep the Slater 
determinants with an energy less than 5 MeV and with 
a width smaller than 1.7 MeV. 

The results of calculations for the 0"*" states in ^^O are 
displayed in Fig. 8. In this case, one obtains two bound 
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FIG. 8. Complex energies of the 0"*" states in re- 
sulting from the diagonalization of the GSM Hamiltonian. 
One- (In) and two-ncutron (2n) emission thresholds arc indi- 
cated. The physical bound and narrow resonance states are 
marked by squares. The remaining eigenstates represent the 
non-resonant continuum. 



states. An eigenstate lying just close to the real energy 
axis just above the two-neutron threshold is a candidate 
for a resonance. Even though the width of this state is 
very small, only the overlap with the states calculated in 
the pole approximation can give the answer. Figure 9 
shows the overlap of the ground-state wave function ^'o of 
^^O, calculated in the pole approximation, with all the 
0+ eigenstates of the GSM Hamiltonian resulting from 
the full diagonalization. One can see that only one state 
(the GSM ground state) has a significant overlap with 
5*0; hence, the identification of the ground state wave 
function is unambiguous. Figure 10 illustrates a more 
challenging case of the third O3" state of ^^O. In spite of 
the fact that this state is embedded in the non-resonant 
0+ continuum, its identification is straightforward. It 
is interesting to notice that those GSM states in Fig. 



9 




Energy of I (MeV) 

FIG. 9. Absolute value of the overlap between the ground 

state wave function v|/o of calculated without the couphng 
to the non-rcsonant continuum (pole approximation) and the 
different 0^ eigenstates of the GSM Hamiltonian (calcu- 
lated with the couplings to the (^3/2 non-resonant continuum), 
as a function of the total energy. 
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FIG. 10. Similar to Fig. 9 except for the O3" state wave 
function of ^*0, *3- 



B. Completeness of the many-body Gamow basis: 
example of two interacting particles 



8 that represent the non-resonant background tend to 
align along regular trajectories. As discussed in Refs. 
[26,27], the shapes of these trajectories directly reflect 
the geometry of the contour in the complex fc-plane. In 
the t-wo-particle case, this information can be directly 
used to identify the resonance states. 

Figure 1 1 shows the results of calculations for the 0"*" 
states in ^"0. As compared to the ^^O case, the num- 
ber of many-body states is much larger and the regu- 
lar pattern of non-resonant states reflecting the struc- 
ture of the contour is gone (the figure represents the pro- 
jection of four-dimensional trajectories onto two dimen- 
sional space). While the two lowest (bound) states can 
be simply identified by inspection, for the higher-lying 
states it is practically impossible to separate the reso- 
nances from the non-resonant continuum. However, the 
procedure outlined above makes it possible to identify 
unambiguously the many-body resonance states. On the 
other hand, the method proposed in Refs. [26,27] cannot 
be easily applied. 

The many-body resonances should be stable with re- 
spect to small deformations of the contour (as the physi- 
cal solutions should not be dependent on the deformation 
of the basis). This observation offers an independent cri- 
terion for identifying resonance states. Figure 12 shows 
the effect of a small deformation of the contour on the 
stability of selected 0+ states in ^"^O. As expected, only 
the states which have previously been identified as res- 
onances are stable with respect to small changes of the 
contour; the states belonging to the non-resonant contin- 
uum 'walk' in the complex energy plane following con- 
tour's motion. 



In this section, we shall discuss the completeness of the 

many-body basis spanned by Gamow states. Since the 
number of configurations is growing extremely fast with 
the number of valence nucleons, and this is enhanced by 
including the non-resonant continuum, we shall restrict 
our discussion to the case of two valence particles. By 
analyzing the behavior of the wave function as the num- 
ber of basis states grows, one can assess the impact of 
various truncations in the valence space. This is espe- 
cially important for calculations with zero-range forces, 
such as the SDl interaction employed in this work, which 
require an energy cut-off. In the calculations contained 
in this section, the configuration space consists of all the 
Slater determinants of energy less than 35 MeV. 

As a first example, we shall consider the convergence 
of the ground-state energy of with the increasing size 
of the non-resonant phase space. The GSM s.p. space 
consists of the 0^5/2, lsi/2 orbitals and the 0^3/2 Gamow 
resonance. This discrete basis is supplemented by adding 
successive continua: si/2, ■ • • j in the decreas- 

ing order of their importance. Since 0(i3/2 is a resonance, 
a (^3/2 contour should be complex, and we take it ac- 
cording to Sec. VA. Other non-resonant continua are 
real and for their path we choose a straight segment with 
< /c < 1.3 fm~^ . The discretization of each continuum, 
whether real or complex, is made with 10 points. 

As one can see in Fig. 13 (top), the convergence is 
achieved with the ^3/2, /7/2, (^5/2, si/2, and j^ji con- 
tours; the contributions from all the remaining partial 
waves with I > 2> and £ = 1 are practically negligible. 
As a rule of thumb, the continua which contribute most 
to the binding energy are those which are associated ei- 
ther with s.p. resonances in the basis or with weakly 
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FIG. 11. Similar to Fig. 8 except for the 0+ states in ^°0. 
The small negative widths of several unbound states are due 
to the fact that the number of discretization points along the 
contour is fairly small (5) in this case. 



bound s.p. states. This is the case with the 0^5/2, lsi/2j 
and 0(i3/2 orbitals, and also with the O/7/2 state which 
has an energy of +7.43 MeV and a width of 3.04 MeV. 
Even though this wide resonance is absent in the basis, 
it has nevertheless an indirect influence through the /7 /2 
contour which 'keeps memory' of its presence. One may 
notice that the presence of the (complex) ^3/2 contour 
affects not only the real part of the ground-state energy 
but also its imaginary part; the ground-state energy has 
a spurious width of ~ -130 keV in the pole approxima- 
tion, and this width is reduced to ~2 keV when the c?3/2 
(10-point) contour is added. The spurious width of the 
ground state remains stable when other (real) contours 
are added. 




Re(E) (MeV) 



FIG. 12. The effect of small changes in the contour on the 
stability of resonant and non-resonant 0^ states in ^"O. Top: 
the contour in the complex-fc plane corresponding to the Od3/2 
continuum. The direction of the contour's deformation is in- 
dicated by an arrow. The calculations were performed for 
four contours, each divided into nine segments (i.e., 10 dis- 
cretization points); only the first and maximally deformed 
contours are shown. Bottom: the resulting shifts in positions 
of many-body states corresponding to the complex energy re- 
gion of Fig. 11 marked by a dotted line. It is seen that the 
states identified as resonant are very stable with respect to 
small changes of the contour, while the states representing 
the non-resonant continuum move significantly in the direc- 
tion indicated by an arrow. 



Nucleus 


Configuration 


c2 


180 


0^5/2 
1 

Od^/2 


0.872-hil.l46 xlO"* 
0.044-15.973 xlO"^ 
0.028-16.624 xlO'^ 






0.042+14.709 xlO"^ 






0.015+11.795 xlO"^ 


«He 


OP3/2 

0p?/2 


0.891-10.811 
0.004-10.079 
0.255+10.861 






-0.150+10.029 



TABLE II. Squared amplitudes of different configurations 
in the ground states of two-neutron systems ^**0 and ''He. 
The sum of squared amplitudes of all Slater determinants, in- 
cluding n particles in the non-resonant continuum, is denoted 
by L^i^' . See text for details. 

The structure of the ground-state wave function of 



^^O calculated with the full non-resonant continuum dis- 
cussed in the context of Fig. 13 is given in Table II. 
One can see that the configuration with 2 neutrons in 
the 0(i5/2 shell dominates; the remaining configurations, 
including those with one and two neutrons in the non- 
resonant continuum, contribute with 15% to the wave 
function. Imaginary parts of squared amplitudes are gen- 
erally very small. This is due to the small width of the 
0c?3/2 resonance included in the basis. 

As a second example, we shall investigate the energy 
convergence for the weakly bound ground state of ^He. 
We shall assume that the structure of ^He can be de- 
scribed in the full Op shell with two valence neutrons. 
Since the GSM s.p. space consists of the 6^3/2, Opi/2 
Gamow resonances, the associated P3/2 and Viji non- 
resonant continua should be complex. Here, the P3/2 
contour consists of three straight segments connected at 
points: fco — 0.0 — iO.O, k\ — 0.3 — iO.35, ki — 0.5 — iO, and 
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Number of {£,]) continua 

FIG. 13. Ground-state energy of (top) and ''He 

(bottom) as a function of the number of different {£,j) 
non-resonant continua in the valence space. For 
we employ the complex ^3/2— contour consisting of three 
straight segments connecting the points: fco = 0.0 — iO.O, 
ki = 0.3 - iO.2, k2 = 0.5 - iO, and fcs = 1.3 - iO (all in 
fm-i). 



ks = 0.7 — iO. In the pi/2 non-resonant channel, straight 
lines join the points: fco = 0.0 — iO.O, fci = 0.4 — iO.5, 
k2 = 0.5 — iO, and = 0.7 — iO. For all other contours we 
take a segment [0:0.7] fm"^ of the real fc— axis. The P3/2 
contour is discretized with 20 points. All other contours, 
including the complex pi/2— contour, are discretized with 
10 points. 

As seen in Fig. 13 (bottom), the full energy conver- 
gence is attained with P3/2 and pi/2 contours; other scat- 
tering waves with i > 1 axe negligible. Coupling to the 
non-resonant and pi^2 continua changes not only 

the ground-state energy but also its (spurious) width. 
In the pole approximation, the calculated ground state 
has a huge width of ^ -2 McV, which is reduced to ^ 
-650 keV when the P3/2 contour is added, and it reaches 
r^-lO keV when both P3/2 and pi/2 non-resonant con- 
tinua are included. Contrary to the case of ^^O, the pole 
approximation is totally unreliable, and it does not even 
give a rough approximation to the energy and wave func- 
tion of the ground state of ^He. The main reason is that 
in this case one attempts to describe the bound state 
using the basis which does not contain any bound state 
and, therefore, the non-resonant continuum is essential 
for compensating the resonance contribution. The anal- 
ogous situation, discussed in the context of a one-body 
problem, can be found in Sect. Ill A (Fig. 6). 

The structure of the ground-state wave function of ^He 
including all the (£, j) continua shown in Fig. 13, is shown 
in Table II. One can see that the configuration with two 
neutrons in the 0^3/2 shell dominates, though the imag- 



inary part of the corresponding squared amplitude is al- 
most equal in magnitude to the real part. The amplitude 
of the OPi/2 configuration is small, whereas the contri- 
butions from one and two particles in the non-resonant 
continuum, lI^'^'' and i^'', are almost equally important. 

VI. GSM STUDY OF OXYGEN ISOTOPES 

In this section we shall discuss the GSM results for a 
chain of oxygen isotopes with several valence neutrons 
{Nyai > 2). In the calculations, we assume the core of 
"'^^O. As discussed in Sec. IV, the valence neutrons are 
distributed over the Lsi/2 and 0^5/2 bound shells, the 
0(^3/2 Gamow resonance, and the discretized ^3/2 non- 
resonant continuum. The complex ^3/2 contour corre- 
sponds to two segments defined by the points: fci = 0-|-zO, 
k2 = 0.2 - 1(0.05, and k^ = 0.4 - iO.O. It is discretized 
with 9 points. Consequently, the discretized s.p. GSM 
space (31) consists of a total of 12 subshells on which va- 
lence neutrons arc distributed. Let us reiterate that we 
also introduce the cut-off in the configuration space of 
the GSM. In this study, we shall include all Slater deter- 
minants with an energy (width) smaller than 5 MeV (1.7 
MeV). Moreover, we shall take only Slater determinants 
with, at most, 2 neutrons in the non-resonant contin- 
uum. We have checked that configurations with a larger 
number of neutrons in the non-resonant continuum are 
of minor importance in all nuclei from "'^^O to ^^O, and 
their contribution is < 0.01% in all cases studied. Our 
aim is not to give the precise description of actual nuclei 
(for this, one would need a realistic Hamiltonian and a 
larger configuration space), but rather to illustrate the 
method, its basic ingredients, and underlying features. 

1. Ground states of oxygen isotopes 

According to our calculations (see, e.g.. Table III for 
^^O), the ground-state wave functions of i8-22q g^j.g (Jqjjj. 
inated by the single shell model configuration Od^™' hav- 
ing a weight of 80-95%; the non-resonant continuum 
contribution is relatively small. A one-neutron contin- 
uum provides between 0.1 and 1% of the wave function, 
whereas a two-ncutron continuum yields a contribution 
which varies from ~ 0.01 to ^ 0.1% in these nuclei. This 
means that the main effect of the non-resonant contin- 
uum in these states is to cancel a spurious width induced 
by the 0^3/2 rcsonancx; included in the basis, and their 
influence on the real part of the energy can be neglected. 
It is interesting to compare the squared amplitudes for 
-'^^O in Tables 11 and 111. The values in Table 11 have been 
obtained by taking all the contours from si/2 to /111/2; 
whereas the results in Table III have been calculated with 
only a d^/2 contour. Increasing the non-resonant space 
gives rise to the depletion of the occupation of bound 
shells. 
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Configuration 

j2~ 



Od; 



Od: 



■5/2 

3/2 
'(1) 



0.935+il.836xl0-* 

0.040-il.560xl0"''' 
0.021-14.973x10"^ 
3.913x10-^+14.412x10"^ 

- 1 .002 X 10"'* +i:!.9:i9 X 10"" 



TABLE III. Same as in Table IT but only including the 63/2 
contour. 



The predicted one-neutron separation energies Sn for 
oxygen isotopes are displayed in Fig. 14 and compared 
with the data. The calculations reproduce experimental 
separation energies, including the magnitude of the odd- 
even staggering. The difference between GSM and the 
data are in the worst case ~ 500 keV. 
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FIG. 14. Calculated (crosses connected by a solid line) and 
experimental (circles with error bars) one-neutron separation 
energies in the 18-22 Q 



2. Level schemes of oxygen isotopes 

Shown in Figs. 15 - 17 are calculated and experimental 
spectra of ^'^"^"O. Let us stress that the main purpose of 
this comparison is to demonstrate the internal coherence 
of GSM results and the numerical feasibility of GSM in 
realistic applications. Nevertheless, in spite of the sim- 
plicity of the Hamiltonian, the overall agreement between 
the GSM and experimental spectra is quite reasonable. 

Figure 15 illustrates the case of ^^O. One can see that 
the first three lowest states are relatively well described. 
Moreover, the one-neutron emission threshold is repro- 
duced. In GSM, this threshold is obtained from the dif- 
ference of calculated ground-state energies of and 
^'^O. It is important to note that all states above the 
one-neutron emission threshold are predicted to be reso- 
nances, i.e., their widths are positive, whereas all states 
below this threshold are calculated to have small spurious 
negative widths. This feature is certainly not imposed on 
the model and it proves both an internal consistency of 
the GSM as well as a correct numerical account of the 
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FIG. 15. The GSM level scheme of calculated in the 
full sd space of Gamow states and employing the discretized 
(10 points) 1^3/2 non-rcsonant continuum. The dashed lines 
indicate experimental and calculated one-neutron emission 
thresholds. Experimental data are taken from Ref. [51]. 



completeness for many-body Gamow states. One should 
notice that in the calculations with a discretized contin- 
uum, it is impossible to assure that the bound states have 
a width exactly equal to zero. Instead, one expects that 
this numerical artifact of continuum discretization be- 
comes numerically unimportant with increasing the num- 
ber of points on the contour (see Table 1 of [25]). With 
10-point discretization, the spurious width remains ^ —5 
keV for all bound states of ^*0, and we consider this a 
reasonable precision. Similarly, for resonances lying very 
close to the lowest particle emission threshold, hence hav- 
ing narrow widths, it may happen that their width results 
from the calculations as negative if the many-body com- 
pleteness relation is numerically not well satisfied (see 
Fig. 11 and related discussion). 

The calculated level scheme of ^^O (see Fig. 16) repro- 
duces the main experimental features. The ground state 
is 5/2+, as observed, but the higher-lying states 3/2^ 
and 1/2^ are given in the reversed order. Nevertheless, 
the next four states are rather well reproduced, in spite 
of the inversion of the 5/22" and 3/22" levels. Finally, the 
GSM prediction for ^"O is shown in Fig. 17. Here, the 
overall agreement between calculations and experiment 
is best for all the isotopes studied. 



3. Distribution of valence particles 

Radial features of the density distribution of valence 
nucleons arc always of great theoretical interest. They 
determine the nuclear size, electromagnetic transition 
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FIG. 16. Same as in Fig. 15 except for ^^O. As the num- 
ber of states becomes large above the one-neutron emission 
threshold, only selected resonances are shown. The electro- 
magnetic transitions listed in Table V are indicated by arrows. 
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FIG. 17. Same as in Fig. 15 except for ^^O. 



decaying state the particle distribution cannot be given 
with an unlimited precision. 



rates, polarization charges, and many other nuclear prop- 
erties. To assess the radial extent of the nucleonic 
distribution obtained in the GSM, wc investigate the 
monopole form factor defined through the normalized ra- 
dial Gamow wave functions Ui{r) (26): 

p{r) =Y,ni\u,{r)){u,{r)\, (39) 

i 

where rij = (a+a„) is the GSM occupation coefficient of 
the s.p. Gamow orbital i. By definition, the form factor 
is normalized to the number of valence particles: 

/ p(r)dr = iVvai, (40) 

Jr=0 

where the "prime" sign indicates that the integral is cal- 
culated by means of complex scaling. 

Figure 18 shows the neutron form factors for the 
ground state and the highly excited resonance state of 
-•^^Q. As for all observables associated with bound states, 
the neutron form factor in the ground state of is real 
even though the 0^3/2 resonant state is included in the 
basis. The imaginary part of p{r) is very small; its value 
is less than 3 x 10~^ for all values of r. On the other 
hand, for the 25" resonance, having a width of ~150 keV, 
the form factor is always complex. This feature is not an 
artifact of the discretization procedure. In fact, the imag- 
inary part of the form factor has an interesting physical 
interpretation. Namely, it is associated with the uncer- 
tainty in the determination of the mean value (given by 
the real part) [34,52,30,53,54]. In other words, in the 



4.. Electromagnetic transition probabilities 

In the shell-model calculations with Gamow states, 
only radial matrix elements are treated differently as 
compared to the standard shell model. This means 
that the electromagnetic (EM) transition selection rules 
and the angular momentum and isospin algebra do not 
change. To calculate the EM transitions, one can no 
longer use the long wavelength approximation because of 
the presence of the non resonant continuum. Indeed, for 
the diagonal EM matrix elements {nilijiki\0\nflfjfkf) 
{ki = kf) between the scattering states, the complex 
scaling cannot be carried out (see discussion around 
Eq. (17)). Furthermore, since in the long wavelength ap- 
proximation the EM operators behave like , one has to 
deal with derivatives of delta functions, which is difficult 
to handle. Without the long wavelength approximation, 
however, these matrix elements become finite, because it 
is always possible to carry out a complex scaling with the 
Bessel function of the photon ^^(gr), asq ^0. Moreover, 
as they represent a set of measure zero, the diagonal non- 
resonant EM matrix elements can be put to zero in the 
discretized calculation. As all the other matrix elements 
can be regularized, the EM matrix elements are all well 
defined. 

Tables IV and V display the selected EM transition 
rates in -'^^O and ^^O calculated in GSM with and without 
the contribution from the non- resonant continuum. 
In all calculations, we have taken the effective neutron 
charge (5e„ = 0.5e. While the radial form factor discussed 
in Sec. (VI 3) determines the structure of monopole oper- 
ators, the EM probabilities probe the off-diagonal matrix 



14 




-0.1 - 



0123456789 10 

r(fm) 

FIG. 18. Monopole form factor of valence neutrons (39) 
in the ground state (top) and the 2g" resonance (bottom) of 
^*0. Real and imaginary parts of the density are shown with 
solid and dotted lines, respectively. According to definition, 
Eq. (40), the form factor is expressed in units of fm~^. 



sitions, the experimental values are reproduced within a 
factor of 3 by the GSM, and we find this agreement sat- 
isfactory in light of a very simple Hamiltonian employed. 



Transition 




Without contour 


With contour 


Exp. 




^0+ 




2.75 iO.OOl 


2.68-10.011 


3.32 




^2+ 




1.99+10.015 


1.94-10.003 


1.19 




-^2+ 




1.66+10.007 


1.48+10.006 
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^2+ 




0.04-iO.OOO 


0.04-iO.OOO 


0.14 


2J 


^0+ 




0.19+10.004 


0.18-iO.OOl 


1.3 



TABLE IV. Electromagnetic GSM rates (aU in W.u.) for 
the selected transitions in calculated (i) in the pole ap- 
proximation without the complex-energy contour and (11) with 
the dz/2 contour representing the non-resonant continuum. 
Experimental data [51] are shown in the last column. 



Table V displays transition rates in ^^O, indicated in 
Fig. 16 by arrows. The three lowest transitions are be- 
tween bound states and the GSM calculations with con- 
tour predict a very small imaginary part (less than 0.001 
W.U.). The next four transitions in Table V involve the 
unbound 7/2^, l/2j, and 5/2;^ levels lying above the 
one-neutron threshold, and the corresponding transition 
rates are, in general, complex. The effect is particu- 
larly pronounced for the 1/2 J 5/23^ transition be- 
tween unbound states. As mentioned in Sec. VI 3, the 
imaginary part gives the uncertainty of the average value 
[30,53,54]. In all cases, the real part of the matrix ele- 
ment is slightly influenced by the interference with the 
non-resonant background. 



elements; hence the continuum coupling manifests itself 
differently. As seen in Tables IV and V, the non-resonant 
continuum plays an important role both for transitions 
involving bound states only and for transitions involving 
an unbound state or states. In all cases, the real part 
of the transition rate is reasonably approximated in the 
pole approximation. 

Similarly, as in the problem of spurious negative widths 
of bound states, the imaginary part of EM transitions 
between bound states must disappear in the limit of the 
non-discrctized continuum. Indeed, in our calculations 
the imaginary part of such a transition rate tends to 
zero when the coupling to the non-resonant continuum 
is added. In the cases studied, taking only 10 points on 
the contour decreases the imaginary part by one to 
two orders of magnitude, depending on the transition. 
The real part of the transition probabilities changes by 
up to ~ 10% when the contribution from the contour is 
considered. This is nevertheless an important change if 
one notices that the configurations with one and two neu- 
trons in the non-resonant continuum amount to ~ 0.1% 
of bound wave functions in these nuclei. For most tran- 



Transition 


Without contour 


With contour 


Exp. 


3/2+ 


Ml 


5/2+ 


0.01-10.002 


0.01-10.0 


0.088 


1/2+ 


Ml 


3/2+ 


0.02-10.003 


0.02+iO.O 


0.0093 


1/2+ 


B2 


5/2+ 


3.07+10.010 


3.07+10.003 


0.58 


9/2+ 


E2 


5/2+ 


1.78+10.000 


1.74 i0.006 


<1 


7/2^ 


E2 


5/2+ 


0.12-10.062 


0.13-10.037 




1/23+ 


E2 


5/2+ 


4.57+i0.430 


4.58+10.373 




7/23+ 


E2 


3/2+ 


0.12+10.038 


0.13+10.054 




l/2j 


E2 


3/2+ 


0.82-0.034 


0.85-10.012 





TABLE V. Same as in Table IV except for electromagnetic 
transitions in '^^O shown in Fig. 16. 



VII. GSM STUDY OF HELIUM ISOTOPES 

A description of the neutron-rich helium isotopes, in- 
cluding Borromean nuclei ^'^He, is a challenge for the 
GSM. ^He is a well-bound system with the one-neutron 
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point No. 



P3/2 



Pl/2 



Nucleus 



0.01-iO.Ol 

0.10-i0.06 
0.20-i0.095 
0.30-i0.06 
0.40-iO.Ol 



0.05-i0.05 

0.20-i0.43 
0.40-i0.45 
0.60-i0.43 
0.75-i0.05 



Configuration 





TABLE VI. Discretized pz/i and pi/2 contours in the 
complox-fc plane (in fm~^) representing the non-resonant con- 
tinuum in the He calculation. 



emission threshold at 20.58 MeV. On the contrary, as dis- 
cussed in Sec. IV, the nucleus '"^Hc is a broad resonance. 
The two-neutron system ®He, on the contrary, is bound 
with the two-neutron emission threshold at 0.98 MeV and 
one- neutron emission threshold at 1.87 MeV. The first ex- 
cited state 2+ at 1.8 MeV in ^He is neutron unstable with 
a width r=113keV. 

The s.p. configuration space used includes both reso- 
nances 0p3/2, Opi/2 and the two associated complex con- 
tinua and ]3i/2) which arc discretized with 5 points 
each. The fact that the resonances in the basis are so 
broad requires particular attention when selecting points 
along the contour. They are given in Table VI. In 
this case, we cannot significantly increase the density 
of points along the contour because the calculation for 
the heavier helium isotopes with more valence neutrons 
would not be feasible. On the other hand, we do not in- 
troduce any restriction on energies and widths of Slater 
determinants included. Also, no restriction on the num- 
ber of neutrons in the non-resonant continuum is im- 
posed, with the only exception being ^He, where we al- 
low for at most 4 neutrons to occupy "contour shells". 
Contrary to the discussion in Sect. VB, we neglect all 
remaining real continua in the present calculation. Wc 
shall see that this will have an impact on the relative 
weight of different configurations in the wave function. 



5. Spectra of helium isotopes 

The energies of the lowest GSM states of helium iso- 
topes are shown In Fig. 19, and the structure of their 
ground state wave functions is given in Table VII. As 
seen in Table VII, the non-resonant continuum contribu- 
tions L^-* are always essential, and, in some cases (e.g., 
^'^He), they dominate the structure of the ground-state 
wave function. Moreover, as can be seen in the example 
of ^He, configurations with many neutrons in the non- 
resonant continuum are essential for fulfilling the com- 
pleteness relation. In this particular case, the contri- 
bution is even more important than the contribution from 
the resonant states, and even the L^'' configuration gives 
a non-negligible contribution of the order of 1-2%. (i^^ 



"He 



OP3/2 
OP3/2 



r(2) 



0.870-i0.736 
0.007-10.076 
0.271 -l-iO.752 
-0.147-FiO.060 



^He 



*He 



OP3/2 
OP3/2OP1/2 

0}33/20}3?/2 

r(3) 

-^f* 

OP3/2 

0P3/2 0P?/2 

L« 

(2) 



l.llO-iO.879 
0.006-10.029 
0.022-10.042 
0.050 +10.951 
-0.185-M0.008 
-0.002-10.009 



L 



0.296-11.323 
-0.060-10.158 

1.596-hil.066 
-0.728+10.630 
-0.125-10.204 

0.020-10.012 



^He 



OP3/2OP1/2 

r(4) 



0.180-11.328 
1.596+10.734 
-0.584+10.801 
-0.217-10.193 
0.025-10.034 



TABLE VII. Squared amplitudes of different configura- 
tions in the ground states of ^~®Hc. The s.p. space consists of 
both 0p3/2, Opi/2 Gamow resonances and the two associated 
complex continua P3/2 and pi/2 which are discretized with 5 
points each. 



amplitudes seem to decrease with n.) Without the con- 
tour, the predicted ground-state energy of ^He is +2.08 
MeV and the spurious width is huge, F = —4.16 MeV. 
The inclusion of scattering states lowers the binding en- 
ergy to —1.6 MeV. Clearly, the pole approximation fails 
miserably in this case. 

It is instructive to compare the wave function decom- 
position in Tables VII (only and pi/2 contours are 
included) and II (all contours up to fj/2 are included). As 
expected, the spread of the ground-state wave function 
over non-resonant continuum states is larger in the latter 
case, but this effect is relatively less important than in 
the case of ^^O (cf. Fig. 13). 

Our calculations reproduce the most important fea- 
ture of ®He and ®He: the ground state is particle-bound, 
despite the fact that all the basis states lie in the contin- 
uum. In spite of a very crude Hamiltonian, rather limited 
configuration space, etc., the calculated ground state en- 
ergies shown in Fig. 19 reproduce surprisingly well the 
experimental data. The odd- A'' isotopes of ^'^He are cal- 
culated to be neutron resonances. The neutron separa- 
tion energy anomaly (i.e., the increase of one- neutron 
separation energy when going from ^He to ^He) is repro- 
duced. One can see that ^He is predicted to be better 
bound than ^He, though, contrary to the data, the two- 
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s.p. basis: ^He 

pj,2 0.75 -iO.3 MeV u„b„„„<j, Vrps = SDI-force 



0* 



Helium isotopes 



*He 'He 



0+ 



8 



He 



^He 



FIG. 19. Experimental (EXP) and predicted (GSM) bind- 
ing energies of ®~^He as well as energies of 3^=2^ states in 
^'*He. The resonance widths are indicated by shading. The 
energies are given with respect to the core of ^He. 



neutron separation energy of ^He is smaller than that of 
^He. Also, the energies of excited 2^ states are in fair 
agreement with the data. 

The present calculations for helium isotopes have a sys- 
tematic tendency to underbind. The spurious width of 
the ground state of ^He is ~ 100 keV, largely due to the 
Op3 /2 resonance. For that reason, the removal of the spu- 
rious width (which reflects the accuracy of our calcula- 
tions) cannot be done here as precisely as in the oxygen 
case. In order to check the stability of the results for 
^'^He, we carried out additional calculations increasing 
the number of points along the complex contour. We 
have found that the value of the two-neutron separation 
energy difference, 5'2n(*He) - S'2n(®He), which is nega- 
tive with 5 states in each non-resonant continuum, de- 
pends sensitively on the number of discretization points. 
This is because S'2n(*He) increases fast with the size of 
the non-resonant phase space. While with our schematic 
Hamiltonian we were unable to find a fully satisfactory 
result for the ground-state energy of ^He just by increas- 
ing the number of discretization points, we can see that 
there is a systematic improvement when the number of 
non-resonant shells increases. Therefore, one would be 
tempted to associate the so-called 'helium anomaly' in 
the position of two-neutron emission threshold when go- 
ing from ^He to ®He with the strong enhancement in the 
occupation of non-resonant continuum states. 



6. Distribution of valence neutrons in ^ He 

Form factors of valence neutrons (39) in 0^- and 2+ 
states of ®He are shown in Figs. 20. In this calculation, 
we have taken discretized pi/2 and continua with 20 
points each. Even with this very high number of contin- 
uum states, the imaginary part of the form factor in the 
ground state of ^He is still ~ 2x 10"'^ fm"-'^. It is inter- 
esting to see that the radial distribution of neutrons in 
this case extends well above the range of the one-body 
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FIG. 20. Monopole form factor of valence neutrons for the 
ground state (top) and for the lowest 2f resonance (bot- 
tom). of '^He. 



potential, as expected for a weakly bound system. In 
the case of the first excited (resonance) 2'^ state, p{r) is 
complex; real and imaginary parts of the density are of a 
comparable size, which is obviously related to the large 
width (~ 500 keV) of this state. 



VIII. SUMMARY 

This paper contains the detailed description of the 
Gamow Shell Model, which can be viewed as a straight- 
forward extension of the standard diagonalization shell 
model that allows for a consistent treatment of bound 
states and the particle continuum, including both reso- 
nances and the non-resonant background. Our first appli- 
cation, based on a realistic, albeit simple, two-body SDI 
force, concerns many-neutron systems, including states 
and nuclei near the neutron-emission threshold. 

In this work, we succeeded in overcoming several obsta- 
cles which traditionally plagued previous continuum shell 
model applications. In addition to the successful inclu- 
sion of the continuum-continuum coupling by means of 
the complex rotation technique (exterior complex scal- 
ing), we incorporated the non-resonant part of the con- 
tinuum. This has been achieved by discretizing the con- 
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tour in the complex fc-plane for each partial wave. An- 
other problem which has been solved in our study is the 
isolation of resonance states. As a result of the GSM 
diagonalization, one obtains a multitude of states corre- 
sponding to the many-body continuum, some being res- 
onances and some representing the non-resonant back- 
ground. Our work offers a simple prescription on how 
to identify the resonance states from the multitude of 
complex-energy cigcnstatcs of the GSM Hamiltonian. 

The GSM Hamiltonian is composed of a one-body term 
and a two-body residual interaction which is directly 
written in terms of space, spin, and isospin coordinates. 
The two-body matrix elements have to be determined 
separately for each case by means of radial integration. 
As a result, the resulting Hamiltonian matrix fully takes 
into account the continuum coupling, in particular the 
spatial extension of s.p. wave functions, determining the 
physics of halo nuclei. 

We demonstrated that the contribution from the non- 
resonant continuum is essential, especially for unbound 
and near-threshold states. In some cases (e.g., *'^He) 
non-resonant continuum components dominate the struc- 
ture of the g.s. wave function. The so-called 'pole ap- 
proximation' (resonant state expansion) breaks down in 
such cases. In addition, the inclusion of the non-resonant 
continuum also impacts results for bound states. 

The results of our first calculations for binding ener- 
gies, spectra, electromagnetic matrix elements, and nu- 
cleonic distributions are very encouraging. In particu- 
lar, pairing correlations due to the continuum-continuum 
scattering can bind the ground states of '''^He with a com- 
pletely unbound basis provided by the s.p. resonances of 
^He. The 'Helium anomaly' (an increase in one- and 
two-neutron separation energy when going from ^He to 
^He) is explained in terms of the neutron scattering to 
the non-resonant background. In all cases considered, 
our calculations yield neutron resonances above the cal- 
culated neutron threshold - a property that is not guar- 
anteed a priori by the formalism. 

Other applications of GSM, including the case of open 
proton and neutron shells, also employing more realistic 
effective interactions, are in progress. We are also work- 
ing on the optimization of the non-resonant part of calcu- 
lations (choice of the contours, distribution of discretiza- 
tion points, etc.). We are convinced that the Gamow 
Shell Model will become a very useful theoretical tool 
unifying structure and reaction aspects of weakly bound 
nuclei. 
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